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I. INTRODUCTION 



It is well known that the four- dimensional Schwarzschild black hole in Einstein grav- 
ity is stable classically under linearized perturbations. Recently, Ishibashi and Kodama [3] 
have shown that this stable behavior extends to hold for higher dimensional cases. How- 
ever, some black strings or branes, which have hypercylindrical horizons instead of compact 
hyperspherical ones, are found to be unstable as the compactification scale of extended direc- 
tions becomes larger than the order of the horizon radius - the so-called Gregory-Laflamme 
instability . The simplest black string would be the five-dimensional Schwarzschild black 
string in Einstein gravity that is a foliation of four- dimensional Schwarzschild black holes 
along the fifth direction. 

Gregory and Laflamme [3] also considered a class of magnetically charged black p-brane 
solutions for a stringy action containing the NS5-brane of the type II supergravity. For 
horizons with infinite extent, they have shown that the instability persists to appear but 
decreases as the charge increases to the extremal value. On the other hand, branes with 
extremal charge turned out to be stable 0]. Since their discovery of such linearized in- 
stability, black strings or branes have been believed to be generically unstable classically 
under small perturbations except for the cases of extremal or suitably compactified ones, 
and the Gregory-Laflamme instability has been used to understand physical behaviors of 
various systems involving black brane configurations as in string theory. Recently, however, 
Hirayama and Kang (Si] analyzed the stability of three types of black string backgrounds in 
five-dimensional AdS space. With or without the presence of a 3-brane, the geometry of 
these black strings in consideration is warped in the fifth direction, resulting in no transla- 
tional symmetry along the horizon. They showed that the AdS4-Schwarzschild black string 
becomes stable as the horizon radius is larger than the order of the AdS4 radius whereas 
Schwarzschild and dS4-Schwarzschild black strings are unstable as usual. It is possible to 
have stationary black string or brane solutions even in four dimensions when a negative cos- 
mo logical constant is present. Interestingly it seems that all known stationary black branes 
in four dimensions are stable. In particular, the case of BTZ black strings has been checked 
explicitly ^J. 

In the context of string theory black branes that Gregory and Laflamme considered are 
those having magnetic charges with respect to Neveu-Schwarz gauge fields only j^. Having 
found some black string systems in which the Gregory-Laflamme instability is absent as 
mentioned above, Hirayama, Kang, and Lee have also considered a wider class of black 
brane solutions for string gravity in order to see whether or not the stability behavior 
drastically changes. Indeed it turns out that the stability of black branes behaves very 
differently depending on the parameter a that specifies the strength of coupling between 
dilaton and gauge fields in the theory as in Eq. (0) below. That is, for magnetically charged 
static black brane solutions in theories of this form j^, there exists a critical value of the 
coupling parameter acr{D,p) to be determined by the full spacetime dimension D and the 
dimension of the spatial worldvolume p of those black p-branes. The case that Gregory and 
Laflamme studied is precisely when a = Qct- Black branes with horizons of infinite extent 
in this case are always unstable as explained above, and magnetically charged NS5-branes 
of the type II supergravity belong to this class. When < a < Ocr, black branes with small 
charge are unstable as usual. As the charge increases, however, the instability decreases and 
eventually disappears at a certain critical value of the charge density which could be even 
far from the extremal point. Magnetically charged black DO, Fl, Dl, D2, D4 branes of the 
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type II string theory belong to this class for instance. When a > aa, on the other hand, 
the instability persists all the way down to the extremal point. Magnetically charged black 
D5 and D6 branes are in this case for example. However it is shown that all black branes 
mentioned above are stable at the extremal point, which might be expected due to the BPS 
nature of extremal solutions in string theory. 

Recently, Gubser and Mitra gave a conjecture about when the classical instability 
of a black brane sets on in terms of thermodynamic stability. This Gubser-Mitra (GM) 
conjecture is a sort of refinement of the entropy comparison argument given by Gregory and 
Laflamme states that a black brane with a non-compact translational symmetry 

is classically stable if and only if it is locally t her mo dynamically stable. As argued by 
Reall ^^l; when we expand classical perturbations in terms of Fourier modes along the 
horizon having translational symmetry, the set of linearized perturbation equations for a 
black brane becomes the Lichnerowitz equation with an additional mass term for the black 
hole on horizon cross sections. Now it can be seen that the existence of a threshold mode 
for instability of a black brane is related to the presence of a negative eigen mode in the 
Euclidean path integral for the black hole system. Consequently, the partition function gets 
an imaginary contribution, implying a thermodynamic instability of the black hole system 
on the horizon cross sections, and vice versa. This interesting relationship between classical 
dynamical and local thermodynamic stabilities has been checked explicitly for various black 
string or brane systems 0;BE3;E;^^- When the translational symmetry along the horizon 
is broken, one can see some disagreements for on set points for instability as shown in the 
stabihty analysis for AdS4-Schwarzschild black strings in AdS space j^. It also should be 
pointed out that this conjecture simply gives the information about when a black string or 
brane becomes stable or unstable. It does not explain or predict other details of classical 
stabihty behaviors |8i]. 

In the present paper, we analyze the classical stability of charged static black brane so- 
lutions for the theory in Eq. in the case that n = D/2 [i.e., a = 0) with a self-dual 
n-form field strength. This case includes black D3-branes in the type II supergravity, and 
is of interest for several reasons. Firstly, note that the geometry of the spatial worldvolume 
of black brane backgrounds whose linearized stabilities have been analyzed so far in the 
literature is flat. As can be seen in Eq. (0]) below, however, black brane backgrounds to be 
considered in this paper do not have fiat spatial worldvolume, but have a warping factor 
multiplied. Such overall factor can not be removed by finding a suitable conformally equiv- 



alent theory as in Refs. [8|, UJj since the background dilaton field is constant in this case. 
Therefore a sort of Kaluza-Klein reduction of perturbation equations for these black branes 
does not give the standard form of the Lichnerowitz equation with an additional mass term 
for black holes on the horizon cross sections as usual. Secondly, in the s-wave perturbation 
analyses in Refs. fluctuations of the field strength for unstable modes could be set to 

be zero consistently partly because black branes are charged magnetically only. However, 
black brane backgrounds in the consideration are charged electrically as well as magnetically 
since the five-form field strength is self-dual. Subsequently, it is not consistent to set s-wave 
fluctuations of the field strength being frozen as shall be shown below explicitly. Finally, 
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the case of black D3-brane is not included in the proof of the GM conjecture by Reall 
Although he suggests some generalization of the argument in Ref. [ll[ could include the case, 
such generalization seems to be non-trivial for the reasons mentioned above. Moreover, the 
covariant action including a self-dual field strength is not known, and so it is not clear how 
the self-duality condition could be incorporated in the proof. Therefore it is interesting to 



3 



see not only how different features mentioned above affect tlie detailed stability analysis for 
black branes to be considered in this paper, but also whether or not the GM conjecture still 
holds for these "dyonic" branes. 

Actually the classical stability of the black D3-brane has been studied by using a notion 
of universality classes recently by Gubser and Ozakin jl^. For s-wave fluctuations they di- 
mensionally reduce the ten-dimensional type IIB supergravity action to a three-dimensional 
one that contains gravity and two scalar fields only, and perform there the stability analysis 
for static perturbations in a certain specific form. However it is not clear whether or not all 
relevant perturbations in the original theory have been covered in such analysis. In particu- 
lar, the field strength is assumed to be frozen in their dimensional reduction. Although such 
approximation would be of little significance for small charge, it might change the stability 
behavior significantly as the charge increases. Finally, we would like to point out that in 
string theory the system of D3-branes is understood best in the context of the AdS/CFT 
correspondence. Hence the details of stability behavior in gravity side might be very useful 
for understanding corresponding behaviors in the CFT side. 

In sectionini we briefly summarize the local thermodynamic behaviors of black D3-branes 
in order to get a hint at the classical stability predicted by the Gubser-Mitra conjecture. 
In section IIIIl we perform the perturbation analysis explicitly and give numerical results. 
Finally, some possible physical implications of our results are discussed. 



II. THERMODYNAMIC BEHAVIOR 



Let us consider the action given by 



d XV— ^ 

d^x 
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e-^(i? + 4(90)^)-— F, 



(1) 



The second form is written in Einstein frame with a = {D — 2n) / \J'2{D — 2). It is known 
that there is no covariant action for low energy type IIB supergravity due to the self-duality 
condition for the field strength F5 = *F5. However, the action in Eq. (0) with D = 10 and 
n = 5 is quite close to the type IIB supergravity with all other form fields set to be zero. 
When the rank of the field strength is n = D/2, one has a = and the field equations for 
the action in Eq. ([Q) are given by 



1 



M 

Rmn 



0, 



1 



^MPi -Pr,-i^N ~ 7~T-^ 9mn- 



4n!" 



(2) 



In the following we consider a gravity theory in which the dynamics of fields is governed by 
the field equations given in Eq. Q with an additional constraint of self-duality for the field 
strength on solutions such as 

F = *F. (3) 

Static black p-brane solutions for this theory are Q 

ds^ = -Udt" + -^ + R^dnl + W6^jdz'dz^ 
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= W (-Ndt^ + Sijdz'dz^) + W-^ (N-^dr"^ + r'^dnl 
0, F = A (e„ + *€„) , (4) 



where 

-2/(n-l) 



N 



1 W=il + rsinhV , U = WN, R^ = r^W-\ (5) 



Here n = D /2 = p + 2, the p-dimensional spatial worldvolume directions are denoted by 
with i = 1,2, - ■ ■ ,p, and e„ is the volume form on the unit sphere dQ'^. The mass and both 
electric and magnetic charge densities are 



M = /c (n + 4 sinh V) , X = \j^^^-^ksmh2fi, (6) 

respectively.^ Notice that the spatial worldvolume of this brane is not flat except for the 
uncharged case {i.e., W = 1). The maximally charged extremal limit is A; ^ and fi oo 
with the mass density (or, fce^^) fixed. 

Now let us consider a finite segment of this black p-brane with unit worldvolume. Being 
regarded as a thermal system, it has entropy and temperature given by 

4 fl — \ 4 

S ~ (cosh yu) r" and T = (cosh /i)^"^ , (7) 

respectively. Here = fc^/^""^) is the horizon radius. The specific heat capacity is given by 

Here one finds that, as /x increases, the heat capacity changes its sign from negativeness to 
positiveness at a certain critical value of /i given by 

sinhVcr = l/2. (9) 

Thus, if the Gubser-Mitra conjecture holds for this system as well, we expect that these 
black p-brane backgrounds become stable classically under small perturbations for /i > 
Defining an extremality parameter as 

A 2sinh2/i 

^ = 1^ = ^.-1.2 (10) 

Amax n + 4 smh /i 



with Amax = Y ('^ ~ one notice that this critical value corresponds to 

= ^2- ^''^ 
For the black D3-brane {i.e., n = 5) this value is q^r = 2v/3/7 ~ 0.495. 



Note here that the value of the charge density is corrected from the one given in Ref. ^ {i.e., X = 
fcsinh2/x/V2). 
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III. PERTURBATION ANALYSIS 



In this section we perform the classical stability analysis for small perturbations of fields 
at the linear level. Under — > + 5(j), F — >• F + 5F, and qmn Qmn + ^MAf, we have from 
Eq. (j2I) linearized perturbation equations given by 

VH<P - V(mV;v)0/^^^ - V(M0V^) (/.^^^ - ^g'^^'h) = 0, (12) 
-{n - l)F^^«[^-^"-VM/i5^ - Vm(/^J' - = 0, (13) 



+ ^-^(25F(*^l^--^"-lF^)P,...p„_, - (n- 



p 



(14) 



2n! 

with the perturbed self-duality condition 

5(F-T) = 0. (15) 

It is very important to notice that the self-duality condition imposed on the perturbed field 
strength in Eq. is not (5F = 15F, but 

1 1 p p 

5Fm^-m„ = {*SF)mi-m„ + -hFM^-Mr, - 7 TVx^^^^Q^ "~^epPi-p„-iMi- m„- (16) 

In addition the Bianchi identity requires 

d6F = 0. (17) 

In order to see whether or not the black p-brane backgrounds in Eq. (ji} are stable under 
small perturbations, we need to check if there is any solution for linearized equations given 
above that is growing in time while it is regular spatially outside the event horizon. In case 
that the black brane is stable, it would be very difficult to show there exists no such unstable 
solution in general. When the black p-brane background is unstable, however, it is rather 
easy to find a certain type of unstable solutions. Notice first that, since the background 
dilaton field is constant, the fiuctuation of the dilaton field can be seen from Eq. ()12|) to 
be completely decoupled, i.e., V^^0 = 0. Hence one can set 6(f) = 0. Since the background 
fields in Eq. are stationary and invariant under translations in directions of the spatial 
worldvolume, one can also assume that 

5F{x>', z') = e^*+™'"'5F(r, 9"^), hMNix", z') = e^'^'""'"' HMN{r, 9"") (18) 

for unstable mode solutions. Here f2 > and {9"^} are angular coordinates for the n-sphere. 
We denote the coordinate system by {x^} = {x'^, z^} = {t, r, 6'™, z"^}. 

For s-wave perturbations that are spherically symmetric in submanifolds perpendicular 
to the p-dimensional spatial worldvolume, one can easily see that the Bianchi identity in 
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Eq. ()17p restricts the form of SF further, yielding the only non- vanishing components given 
by 



5F,,,i...,p=e^*+^'"^^'g(r). (19) 

In order to see the presence of unstable mode solutions, it suffices to consider the threshold 
modes {i.e., Q = 0) only as usual jl, [l^. Note that we have not imposed any gauge 
such as the transverse gauge {e.g., V Mh^^ = V^/i/2) for metric perturbations Hmn in 
Eq. (I14j) so far. For unstable threshold modes one may set Ht^ = Hu = by choosing a 
suitable gauge as in Refs. [l^ , instead of the usual transverse gauge as in Gregory and 
Laflamme [3|. Now the remaining non- vanishing components of Hmn are Hu, Hrr, Hri, Hgg 
and Hij. By looking at the form of perturbation equations for Hij and Hri, one may have 

= ^C(r) + 6]pir) (20) 



and 



rn • 

Hl = i^r^{r), (21) 



where = X)f=i^i- If P = 0, this is exactly the same ansatz used by Gregory and 
Laflamme j^. The presence of the new term p above can be understood as follows. Under 
the diffeomorphisms associated with a vector field = e^"^^^' C.'^'^r) , we have 

fW'\ 

Hij = Hij - imi^j{r) - imj^i{r) - SijV I — \ ^^(r). (22) 

Since W = drW does not vanish for charged black branes, the term p should appear in 
general. On the other hand, the H^ component transforms as 

W 

H„ = H„ - zmMr) - ^r) + -^Ur). (23) 

Therefore, when the spatial worldvolume of a black brane is a warped fiat space as in the 
present case, one sees that only two of three functions (, p and t] can be eliminated by 
choosing suitable functions ^r{r) and ^(r) in ^j(r) = imi^{r). We set C = = in our 
analysis. Consequently, for threshold modes the non-vanishing components of s-wave metric 
perturbations finally become 

= e'^^''di^g(Hl,Hl,HlZ,Hi) 

= e^™'"'diag {ip,^,x,---,X,P,---,p) , (24) 

where all unknown functions (p, ip, x, and p are functions of the radial coordinate r only. 

In Appendix A, all linearized perturbation equations in Eqs. (|13p and (jl4p are shown 
in components. We have seven coupled ordinary differential equations for seven unknown 
functions ip, ip, Xi C; Vi P and Q. Although we have not used the transverse gauge condition, 
note that in the Reall gauge {i.e., C = V = ^) linearized perturbation equation Eq. ()A6|) 
becomes equivalent to the 2;*-component of the transverse gauge Eq. ()A10|) . As can be seen in 
Eq. ()A11|) . however, Eq. ()A5|) differs from the r-component of the transverse gauge Eq. ()A9|) 
as the black brane gets charged {i.e., W ^ 0). This property is different from the cases 
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studied in Refs. j8|, |ll| where some of linearized perturbation equations in the Reall gauge 
become precisely equal to the transverse gauge. 

On the other hand, the only non-trivial component of the self-duality constraint in Eq. ()16|) 

is 



2Q - X{yD + ij + {n - 2)p - nx + = 0. 



(25) 



It is very interesting to see that this equation is equivalent to the linearized equations for 
gauge fields in Eq. (fT^ as can be seen in Eqs. ()A7|) and ()A8|) . Thus the self-duality constraint 
is automatically satisfied at least for any s-wave threshold mode solutions of the linearized 
equations. Eq. ()25|) shows that the fluctuation of the field strength is proportional to the 
charge density when the brane is charged weakly. It also implies that 5F ~ Q cannot be 
frozen to be zero for a finitely charged brane since the term in parentheses in Eq. (|25|) would 
not vanish in general. For the case of the DS-brane with the gauge ( = rj = for example, it 
is enough to see that the zeroth-order equation in A {i.e., ipo + ipo — 5xo = 0) is inconsistent 
with other uncharged linearized equations. 

In our gauge chosen as above we now have a set of seven equations for five unknown 
functions if, ip, Xi P ci-nd Q, which consists of two zeroth-order, one first-order, and four 
second-order coupled ordinary differential equations. By doing lengthy calculations one 
can show that two of them are not independent of others, resulting in five independent 
equations for five unknowns. From now on we restrict our consideration into the case of 
the black Z)3-brane, i.e., n = D/2 = p + 2 = 5. We expect that the analysis for other 
cases of black branes in Eq. (0)) can be done similarly and that their stability behaviors are 
essentially same as this case. By eliminating the functions P and Q, we end up with two 
second-order semi-decoupled differential equations for ip and ip as follows: 

(r^ - l) if" + mV2(r^ + 6)^(5r^ - l) - SSr^" + (376 + AQy^ + (636" 

-626 - 19)r^2 + 36(56^ - 186 + ll)r^ + 6^(66 + 31)r^ + 36 
r^{r^ - l)(r^ + 6)^ -m2(r^ + 6)^[7r^^ - 2(66 + 5)r^2 - (ll6^ - 126 + 5jr 
+26(36 - 8)r^ - 36^1 + 16r^ Ub + 2)r^2 - 36(36 + 2)r^ - 36(26^ - l)r^ - 6^ 



m 



rK 



7r^ + (6 - 5)r^ + 6 (6 + 2)r^ - 56(6 + l)r^ - 6^ tp' + 



M - l)K 



m 



- »r 



(r^ + 6) [7(6 + 2)r^2 _ (^111^2 + 35 + io)r^ + 6(146 + 9)r^ + 6 
+56 - 2)r^2 + 6(26^ - 136 - 6)r® + 6(26^ + 146 + 3)/ + 6^(26 - 1 



106^ 
^ = (26) 



and 



(r^ - 1) ^" + rr?r'^(r^ + 6)^(5r^ " ^ " + ^^^''^ + ^^^^^^ + (^^^^ 



3586 - 37)ri^ - 6(76^ + 386 - 163)r*^ - 36^(66 + 5)r^ + 6^ 



(r4 - l)K 



m 
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(r^ - l) (r^ + + m^(r^ + fe)^ [Zr^^ - 2(306 + ll)r^^ + (2lb'^ + 1406 + 23 



26(76 + 32)r^ + 6^ 

6(46^ - 86+ 15)r^ - 36^ 
2 



6(46 + - (26^ + 636 + 4)r^^ - h{ 26^ + 136 - b2]r^ 



V^-— r^p-(6 + 3)r^-6 
rK L 



- (116 + 3^ 



+36(26 + 5) 
-6(106 + 9) 
+6^(26 + 3)" 



+ 8r' 



2 4/ 

m r 



(/ + 6) ^ p2 - (116 + 14)r^ + 3(26^ + 46 + 3 



6r^^ + (156 - 4)r'^ - 6(116 + 32)r^ - 6(26 + 3)(46 - 5)r' 



V? = 0. 



Here 



ir(r) = (r^ + 6) [mV2(r4 + _ (^^^ _ 1) (7^^ 



146r^ 



(27) 



(28) 



and the parameter 6 is defined by 

6 = sinh^ fi. (29) 

In writing Eqs. (j26|) and pTjl . we have used that the background metric in Eq. is invariant 
up to an overall multiphcative factor under rescalings of r — >• r/a, — > ru/a, t ^ t/a, 
— >• z^/a for arbitrary constant a. Thus it should be understood that the dimensionless 
quantities in equations above are actually r ^ f = r /th and m ^ rh = r„m with = 
k^/^^~^\ respectively. 

In order to see whether black D3-branes are unstable or not, now we need to check if these 
two equations allow any spatially regular solution outside the horizon for certain values of 
the parameters m and 6. Let us first consider the boundary behavior of the solutions near 
the horizon and at the infinity. Since Eqs. (j^UI) and ()27|) are second-order coupled linear 
differential equations for ip and if), there are four linearly independent mode solutions in 
general. The asymptotic solutions at spatial infinity ( i.e.,r ~ 00) are given by 



^00 

^00 



(r) ~ e±'™r-5/2 ^ + 2^1 + 
V 8m r 
35 1 



+ 



m 



]m? r 



(30) 



up to overall arbitrary multiplicative constants. Note that only half of these mode solutions 
are regular ones. Similarly, in the vicinity of the event horizon [i.e., r = 1 + A), we have two 
regular asymptotic mode solutions and two singular ones. By being regular we mean that 
the solution should not produce any curvature singularity at the horizon j^. Asymptotic 
regular mode solutions are given by 



1 



1 

96 



(6 + 1)W - 16(106 - 3)(6 + lym^ - 64(6 + 1)[6(1096 - 232) + 75]m' 



+81926(26- 3) 



A2 + 



'^ir„ ~ l + ^(6 + l)W (6 +1)W + 32(6-1; 



A + 



(31) 
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and 

^IIr„ ~ A + — 

i^nr„ ~ - — (6+1) 



{h + 1)W - 2(276 - 1)(6 + 1)W - 16[6(196 - 82) + 19] 



A2 + 



1 + 6)^m^ + 8(36 - r 



+16[6(236 - 122) + 23] 



A2 + 



(6 + 1)W + 2(96 - 11)(6 + 1)W 

(32) 



Here 



Li = 192(6 + 1) [(6 + ifm^ - 8(36 - 1) 
Note that these asymptotic solutions become divergent as Li — > 0, i.e. 

~ 8(36-1)7(6 + 1)1 



(33) 



(34) 



To avoid this divergence for such values of parameters m and 6 (or, close to them), one may 
also use another set of asymptotic solutions by linearly combining these solutions. They 
may be given by 



(6 + Ifrrf - 16(106 - 3)(6 + 1)W - 64(6 + 1)[6(1096 - 232) + 75]m' 



+81926(26 - 3) 



A2 + 



Li + 96(6 + l)^m^ (6 + 1)W + 32(6- 1) 



A + 



(35) 



and 



^11 TH 



Here 



1 



1 



(6 + if 771^ - 16(106 - 3)(6 + 1)W - 64(6 + 1)(6(1096 - 232) + 75)m' 



+81926(26 - 3) 



A 



7(6 + if 771^ + 16(86 - 21)(6 + if 771^ - 64(6 + 1)(6(376 + 24) - 77)m' 



3^4 



+81926(26 - 3) 



A 



16 



+ 1)W - 2(6 + 1)^(276 - l)m^ - 16(6(196 - 82) + 19) 



(36) 



(37) 



The first pair is obtained simply by multiplying Li and the second pair by linearly su- 
perposing two solutions such that the term proportional to A^ in is regular, i.e., being 
zero. 

As in Ref. jsf, let us consider two mode solutions whose asymptotic behaviors near the 
horizon are those in Eqs. (jHT|) and (jH^ as follows: 



Vii 



for r 



AjC + Bie"^^v+ for r 

i'liru for r 



' H 
00, 



A,,e "^^v- + i?rre™'"f+ for r ^ 00. 



(38) 
(39) 
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Now any mode solutions that are regular at the horizon can be written by 

^ = C^j + Eifiji, ^ = Ci)i + Eipij. (40) 

At the spatial infinity, they will behave like 

if ~ {CAj + EAjj) e-™^M_(r) + (CBj + EBu) e™^M+(r), (41) 
il) ~ (CA, + EA,^ e-™"M_(r) + {CB^ + EBn) e'^'u+ir). (42) 

In order for these solutions to be regular, coefficients of the exponentially growing parts 
should vanish 

CB, + EBjj = 0, CBr + EBj, = 0. (43) 

The condition that there exist any non-trivial coefficients C and E satisfying Eq. PHj) will 
be 

D{m, b) = BjBj, - BuBj = (44) 

Namely, the existence of unstable mode solutions depends on whether or not there are certain 
values of parameters m and b satisfying Eq. (j44p . 

We have checked this numerically by using Mathematica. In more detail, having given a 
certain form of initial data as in Eqs. (j38p and (j39|) near the horizon, we solve the coupled 
equations Eqs. (j26p and (j27p numerically, and evaluate the determinant D{m,b) by using 
numerical values for ipij^ and ipiji at sufficiently large r. For a given parameter b of black 
D3-brane, we vary Kaluza-Klein mass m only and search for the m* at which D{m*, b) = 0. 
This can be achieved by finding out a m* around which D{m,b) changes its sign. If there 
exists such m*, the corresponding solution is indeed the threshold unstable mode. 

In the method described above, however, it should be pointed out that in addition to 
real solutions there appear some fictitious solutions in parameter space (m, b) due to the 
divergence behavior of initial data in Eqs. (|3ip and (j32p for a certain range of parameters 
as in Eq. (j34j) . One can find that these fictitious solutions exactly coincide with the curve 
defined by Eq. (jnH). By using other set of initial data such as those in Eqs. (fH^ and 
it turned out that such fictitious solutions never appear. 

Behaviors of some threshold masses we have obtained numerically in parameter space are 
illustrated in Fig. ^ and Fig. |21 The diagram on the left hand side in Fig. ^ shows how the 
dimensionless threshold mass m*{b) varies as the "charge" parameter b increases. Note that 
the uncharged black D3-brane is simply a product of the seven-dimensional Schwarzschild 
black hole and three-dimensional fiat space and that our s-wave perturbation analysis be- 
comes equivalent to that of a pure gravitational perturbation for such Schwarzschild black 
3-brane since the gauge perturbation Q{r) becomes zero as the charge density A vanishes. 
In Fig. ^ the dimensionless threshold mass for the uncharged black Z)3-brane turns out to be 
m*(6 = 0) ~ 1.85. We find such numerical value coincides with those obtained in Refs. 0, |^ 
that, in the uncharged limit, can be regarded as the result for the simple Schwarzschild black 
3-brane. 

As the "charge" parameter b increases, the dimensionless threshold mass decreases mono- 
tonically, and approaches zero at 6 ~ 0.50. We have checked numerically that there is no 
solution for larger values of b than this critical value up to about 6 ~ 10 {i.e., q ~ 0.93), 
which is close to the extremal point q = 1. By using the scaling property explained before 
we can obtain the actual threshold mass from our results as 

m*{M,X)=m*{b)/r„= (^^^ m*(r, = l,q), (45) 
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FIG. 1: Behavior of dimensionless threshold masses for black D3-branes at various values of h. 
Note that fh*{h = Q) ~ 1.85 for the uncharged case. On the right hand side the threshold mass for 
the mass density M = 5 is plotted in terms of the extremality parameter q. The extremal point is 
g = 1 and the existence of solutions has been checked up to g ~ 0.93 (i.e., h ~ 10). Note that the 
critical value at which instability disappears is h^t — 0.50 (i.e., gcr — 0.49). 

where m*(r^ = l,g) = {h + Ahy/'^fh* {h) /h^^^ and h = b{q) should be understood as a function 
of M and A through Eq. . The diagram on the right hand side in Fig. Q shows the 
behavior of the threshold mass for the mass density M = 5 fixed as the extremality parameter 
q increases up to g ~ 0.93. Therefore our numerical results show that, although black D3- 
branes with a given mass density are classically unstable under small s-wave perturbations 
when they are charged weakly, as they get charged further this instability decreases down 
to zero at a certain point far from the extremal one. It can be also seen that the critical 
value of q {i.e., gcrNum — 0.49) at which the instability disappears agrees very well with the 
predicted one {i.e., q^r = 2-\/3/7 ^ 0.495) in Eq. (fTTjl through the GM conjecture. 

Now let us consider the extremal black D3-brane, i.e., /c — *■ and /i (or, 6) — > cxd with kb 
(~ ke'^^/A) = M/4: = c fixed. The extremal black D3-brane is expected to be stable since 
they correspond to the BPS ground state in string theory. Our numerical result obtained 
up to g ~ 0.93 {i.e., b ~ 10) above also seems to indicate that there would not appear 
instability mode when we continue our analysis further up to the extremal point q = 1. 
However, 6 ~ 10 is still small compared to 6 ~ oo, and our analysis based on rescaled 
variables is not appropriate to the case of extremal limit. In particular, the initial data 
in Eqs. (j31|) and (j32p cannot be kept small as 6 — oo. Thus we study the extremal case 
separately. By recovering the parameter k and taking the extremal limit in Eqs. (j26j) and 
(|27p. perturbation equations for the extremal D3-brane are given by 



^ P 



5cW + 5r^^( - 7 + mV^) + 50^^(3 + 4mV2) + cr^° (37 + 20m 



2 2 

r 



+3cV6(21 + 10m' 
+cV2fll + 3mV2 



2 2 

r 



, m'^{c + r^Y 



P 



cW + r^°( - 7 + mV^) + 3cr^(4 + 



2 2 

m r 



2cr^ 



he' + 34cr* - 7r' 



2cr^ 



P 



llc^m^ - 7mV^^ 



+cV2 ( 16 + 15mV^ )+ cr'' ( 80 - 3mV 



.2^2 



.2^2 



1p = Q 



(46) 
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m*(M, A) 

1 . 75 
1.5 
1 .25 

1 

. 75 
0.5 
0.25 

2 4 6 8 10 12 14 

M = 5, 10, 15, 20 

FIG. 2: Threshold masses vs. charge density A for black D3-branes having various values of 
mass density M. Extremal values of the charge density are Amax = Mj^pl ~ 3.5, 7.1, 10.6, 14.1, 
and critical values corresponding to 6cr = 1/2 are Acr = \f^Mll ~ 1.8, 3.5, 5.3, 7.0 for M = 
5, 10, 15, 20, respectively. 




and 



^ P 



+cV^(l3 + 30m 



5cW + r^^( - 49 + SmV^) + cV2( - 7 + 20mV2) + cr^° (l87 + 20m 



2 2 

r 



2 2 

r 



^ P 



+ m^r^^{j + mV^) + cWr^ (21 + 5mV^) 



+cr^2 ( _ 384 _ 46m^r^ + Sm^r'' ) + 2cV* ( 16 - 46m^r^ + 5m*r 



.2^2 



.4^4 



,2^2 



,4^4 



+2cV(l6 - 9m + 5m 
2m^r^ 



4 4 

r 



- 6c^ + 17cV^ - 12cr* + r^^ 



P 



+ r^) (6c^ - llcr^ + r^) = 0, 



where 



P{r) = r^(c + r^) c^m^ + ( - 7 + mV^) + cr^ (l4 + 3mV^) + cV^ (l + 3m' 



2 2 

r 



(47) 



(4J 



and m and r are all dimensionful variables. 

Differently from other extremal p-brane cases in Ref. Q , these equations at the extremal 
point are not decoupled. Repeating similar analysis as in the cases of non-extremal black 
Z}3-branes, for given mass density M (or, c) we have checked for various values of m, but 
found no regular solution. It confirms that the extremal D3-brane is stable, at least under 
s-wave perturbations classically. 



IV. CONCLUSIONS 

To conclude, we have investigated the classical stability of black D3-branes in the type 
IIB supergravity under small perturbations. For s-wave perturbations it turns out that black 
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FIG. 3: Threshold masses with respect to the horizon radius for a given charge density A = 2\/2. 
Black dots denote our result and gray dots that of Ref. [l^l- = corresponds to the extremal 
-D3-brane, and the critical horizon radius is r^cr = (2/\/3)^^* — 1.04. 



D3-branes are unstable when they have small charge density. As the charge density increases 
for a given mass density, however, the instability decreases down to zero at a certain value 
of the extremality parameter {i.e., gcr — 0.49), and then black D3-branes become stable all 
the way down to the extremal point. It has also been shown that such critical value at which 
its stability behavior changes agrees very well with the predicted one by the thermodynamic 
stability behavior of the same system through the Gubser-Mitra conjecture. Therefore, 

to the case of black 
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although the generalization of Reall's proof for the GM conjecture 
D3-branes seems to be non-trivial as explained before, our direct comparison confirms that 
the GM conjecture presumably holds to this case as well. 

The peculiar property in this theory that the five-form field strength should be self- 
dual is imposed for linearized perturbations as an additional constraint on the solution 
space. Interestingly, however, it turned out that, for s-wave perturbations, such constraint 
is automatically satisfied since this constraint equation is precisely equal to one of linearized 
equations for the gauge field. It should also be pointed out that the s-wave instability for 
black D3-branes involves a non-vanishing gauge field fluctuation. This property differs from 
other cases of magnetically charged black j9-branes studied in Refs. j3l, \W\ where allowed 
s-wave unstable solutions are for fluctuations of dilaton and gravitational fields only with 
the gauge field frozen. Although our stability analysis described in this paper is restricted 
to s-wave perturbations so far, we expect that the essential stability behavior of black D3- 
branes would not change even if we consider non-s-wave perturbations further. It follows 
because the strongest instabiUty is expected to be carried by s-wave perturbations |3|].^ 

Now let us compare our result with that of Ref. ^] obtained by using the notion of 
universality classes. Since the gauge field is frozen in the dimensionally reduced action in 
Ref. [ij], it is of interest to see how this difference in the stability analysis affects to the 



See also some other argument in Ref. |3 supporting such expectation 
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result. In Fig. |21 by using Eq. ()45|) and the relation between b and the horizon radius given 



the threshold mass in our result is plotted with respect to the horizon radius for a given 
charge density {e.g., X = 2\/2) together with that in Ref. [l^. Here rjj = corresponds to the 
extremal D3-brane, and the critical horizon radius at which the instability sets on {i.e., her = 
1/2) is expected to be r„,, = (2/^3)^/^ ~ 1.04. Since the non- vanishing gauge fluctuation is 
proportional to the charge density for small charge, we expect that the assumption of frozen 
gauge fluctuation in Ref. jl^ niight be flne when the charge density of black D3-branes 
is small compared to the mass density. Thus threshold masses in two methods probably 
agree very well at least for large horizon radii as can be seen in Fig. El Indeed = 1.6 
corresponds to g ~ 0.1 already. Interestingly, however. Fig. El shows that both results are 
still in very good agreement even for horizon radii near the critical value Th ~ 1.04 {i.e., 
q ^ 0.49). Thus we see that the approximation (5F = taken in Ref. does not change 
the results much even for rather "strongly" charged branes. 
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l/(n-l) 
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APPENDIX A: LINEARIZED PERTURBATION EQUATIONS 

For s-wave perturbations in the form of Eqs. (fT^ and ()24|) without imposing ( = rj = 0, 
the hnearized equations in Eqs. and ()14|) can be written in components as follows: 



(tt)-component: 



U^" + -(\nr^^U^W-')'v' + 



2 V 



m 

w 



2 ^ 



U' 



U' 

+ (n-2)-p' + (n-2)^ 



/ AW" 



2^2n 



, [/' , AW" 



U', AW", mf/' Al^"^ ^ 

P+yC' + Y:^C + ^.-Y^Q = o (Ai) 



(rr)-component: 

f/ + nx" + {n- 2)p" + C"] - - [inr^'^UW-^)' ^Ij' - 



2r2" W J 



'AW" m- 
+ 



^ + + 7 [^^r'UW-') nx' + -^nx +{n-2)- [\nUW') p' 



2 ^ ' 2r2" 

AW" f/ 



2^2n 



+("-2)^P+y(tafM'^)'C' + 



2r^ ' 2 
-component: 

Ux" + - (lnr^"f/W-""2)';^' + 



AW" 
2r2" 



C + TTT 2r/' - — r/ - 



AW^ 

»,2n 



■Q = (A2) 



+ - (inrW-i)' (lnr2("-2)f/2)' 



-(n-2) 



f/- 



AW" 



m 

2r2" " "W^ 

rW-i)" . f/ 
2 



rW-i ■ 2 
X + -(lnrW-i) if' 



_U 

2r2" ^ ~ 2" 



(inrW-i)V' 



^ + ^ (inrW-i)' (lnr2("-2)f/2)' + ^ + {n - 2)- (inrW-i)'p' 



^2|y 

AW" 



2^2n 



-(^-2):^p+^(lnrW"i) C 



AW" 
2r2" 



C+^(lnrW-y, + ^g = (A3) 



2" 

A_l^" 



(z* = Z'')-component: 

[/p" + ^(lnr2"[/W"-^)'p' + 
t/iy' , AW" f/iy' 



' W" uw^ 
W ^ '2VV 



(lnr2"f/W-^)'+ (n-4) 



AW" 



2^2n 



m 
W 



H V> + 

2 2r2" 



ly" U W 



+ ^^flnr2"f/W-^)' 

9 ly V y 



AW" 
2r2" 



2 py 



t/iy' , AW" f/iy',, AW", mw \w ^ 
+ 2 Ty + 1^^^ + 2 Ty ^ + + ww"^ ~ = ° 

(r2;*)-component: 

U [if' + nx' + {n- 3)p'] + (inf/Py-^) V - ^ (lnr2"[/iy-^) V 

■^^ + ?^(lnr-f/W-- 



(A4) 



+ — flnrW"^Vny + — 
2 ^ ^ m 



W 2 W 



„2n 



ri = (A5) 
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(z* 7^ z-')-component: 

2r]' + (in r^^'W-^y 7] 



5F: 



' W" UW f o 9 



W ^ 2W 



m 



— [V9 + + nx + (n - 4)p] = 



(A6) 



\[ip + i)+{n-2)p-nx + C]-'^Q = 
A + ^' + {n- 2)p' - nx' + C] - 2Q' = 



Two non-trivial components of the equation V uh^'^^ = |V^/i can be written as 



(r)-component: 



U [if' + nx' + {n- 2)p' + ('] + U'ip~U (inr^^f/l^-^y ^ 
+U (in r'W~')'nx + [(n -2)p + C] + ^V = 



(A7) 
(A8) 



(A9) 



(z*)-component: 



m[ip + + nx + {n - 4)p - (] - 2r( - (inr^"!^"^) 77 = 



(AlO) 



Note that Eqs. ()A5|) and ()A6|) do not contain the function Q{r\ and also that they are 
related to Eqs. (IX9|l and (jAlOjl as 



y(ra + |(iaoD'-^^b-V^ + ^x + (^-2)p] 



2" 



T] - 



mU W 



C 



U 



w 



+ -— (Inr2-^7W 



2m VT' 



U AW" 



n2n 



(All) 



(A12) 



where {Al) with / = 5,6,8,9 denotes the left hand side of Eq. {Al). Thus one finds that 
in the Reall gauge {i.e., C = ^7 = 0) the linearized perturbation equation Eq. ()A6|) be- 
comes equivalent to the 2;*-component of the transverse gauge Eq. ()A10|) . As can be seen in 
Eq. ()A11|) . however, Eq. ()A5|) differs from the r-component of the transverse gauge Eq. ()A9|) 



as the black brane gets charged (z.e., W 7^ 0). This property is different from the cases 
studied in Refs. jl, [ll| where some of linearized perturbation equations in the Reall gauge 
become exactly equal to the transverse gauge. 
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APPENDIX B: PERTURBATION EQUATIONS FOR THE D3-BRANE 

For black DS-brancs (i.e., n = D/2 = p + 2 = 5), the linearized perturbation equations 
above in the Reall gauge ( = r] — become 

(tt) -component: 

1 \ „ 5r^ + 66 + 1 , _ m2(ri2 + 36r« + SfcV^ + b^) - 46(6 + ly 

{b + 2y + b _ Abjb + l) {b + 2y + b 4b{b + l) 

r9(l + Vr^) ^ rio(l + ^ ^ r9(l + b/r') ^ rio(l + b/r'^y^ 

Ab + 2y + b Ab{b+1) _ V2b(6 + 1) ^ 

(rr)-component: 

, , , 5r^ + {4b-3y -2b . m\r^ + bf + 4b{b + iy, 

(6 + 2)r4 + 6 , 46(6 + 1) 2r^ + br^ + b , 46(6 + 1) 

^^r9(l + 6/r4)^ + rio(i + 6/r4)2'^ ^ r9(l + 6/r4) + rio(l + 6/^)2^^ 



(36 + 2)r^- 6 46(6 + 1) 2^26(6 + 1) 



r9(l + 6/r4) ^ rio(l + 6/r4)2 ^ r^^{l + b/r^Y 
(^^) -component: 

1 \ lOr^ + (56 - 6)r4 - & / _ 1 



rV r\l + b/r^) ri2(l + 6/r4) 

-8r^° - 166r^ + 46(36 + 5)r 



m2(r^2 + 36r« + 36V^ + 6^) 



1 - 1/r^ 46(6 + 1) 1-1/r^ , 

^ + r(l + 6/r4)'^ " rio(l + 6/r4)2'^ " r(l + 6/r4)'^ 



2r^ + 46r^ + 6(6 - 1) 1 - 1/r^ , _ 46(6+1) 
rio(l + 6/r4)2 '^+r(l + 6/r4)^ r^^{l + b/r^Y^ 

2^26(6+1) 



Q = (B3) 



rio(l + 6/r4)2 

(2;* — 2;^)-component: 

1\ „ 5/-- + (86- l)/-i -46 , r/;2(,-i2 + 3/,,.s + + //')- 126(6 +1)/' 

~ ^ ^ r9(l + 6/r4) ^ ri2(l + 6/r4)2 ^ 

1-1/r^ , 46(6 + 1) 1-1/r^ , 46(6 + 1) 

r5(l + 6/r4)'^ + rio(l + 6/r4)2'^ ~ r\l + b/r^y ~ rio(l + 6/r4)2^ 

l-l/r'' , 46(6 + 1) 2y^26(6+ 1) 

^^5(l + 6/r4)^^ + rio(i + 6/r4)2^^ " rio(i + 6/^)2^ = (^^) 
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(r2;* )-component : 



, .s 5x = (B5) 



r9(l-l/r4)(l + 6/r4)'"^ r^{l + b/r^) 
(z* 7^ z-')-component: 

<^ + V' + 5x + p = (B6) 

6F: 

A [(^ + ^ + 3p - 5x] - 2g = (B7) 
Here b = sinh^ /z and we set k = 1. 
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